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The calculation of the full electroweak 0(a) corrections to the charged-current four-fermion production pro- 
cesses e + e~ — » v t t + n~Vfj., ud^Tv^, and udsc is briefly reviewed. The calculation is performed using the complex- 
mass scheme for the gauge-boson resonances. The evaluation of the occurring one-loop tensor integrals, which 
include 5- and 6-point functions, requires new techniques. The effects of the complete 0(a) corrections to the total 
cross section and to the production-angle distribution are discussed and compared to predictions based on the 
double-pole approximation, revealing that the latter approximation is not sufficient to fully exploit the potential 
of a future linear collider in an analysis of W-boson pairs at high energies. 



1. INTRODUCTION 

At LEP2, W-pair-mediated four-fermion (4/) 
production was experimentally explored with 
quite high precision (see Ref. £Q and references 
therein). For LEP2 accuracy, it was sufficient 
to include corrections in the so-called double- 
pole approximation (DPA), where only the lead- 
ing term in an expansion about the poles in 
the two W-boson propagators is taken into ac- 
count. Different versions of such a DPA have 
been used in the literature [21314151b') . Although 
several Monte Carlo programs exist that include 
universal corrections, only two event generators, 
YFSWW 03 and RacoonWW (5T7IB| . include 
non-universal corrections. 

In the DPA approach, the W-pair cross sec- 
tion can be predicted within ~ 0.5% (0.7%) in 
the energy range between 180 GeV (170 GeV) and 
~ 500 GeV, which was sufficient for the LEP2 ac- 
curacy of ~ I %. In the threshold region (y/s < 
170 GeV), the DPA is not reliable, and the best 
available prediction resulted from an improved 
Born approximation (IBA) based on leading uni- 
versal corrections only, and thus possesses an in- 
trinsic uncertainty of ~ 2%. 

At a future International c + c~ Linear Collider 



(ILC) |9I10I11| . the accuracy of the cross-section 
measurement will be at the per-mille level, and 



the precision of the W-mass determination is ex- 
pected to be ~ 7MeV from a threshold scan of 
the total W-pair-production cross section |9I1()| . 
For the cross-section prediction at threshold the 
theoretical uncertainty of the IBA is ~ 2% and 
thus definitely insufficient for the planned preci- 
sion measurement of A% in a threshold scan. 

Recently we have completed the first calcula- 
tion of the complete 0(a) corrections (improved 
by higher-order ISR) for e + e~ — > v T r + pT v^, 
ud/i - ^, and udsc, which are relevant for W-pair 
production. We have presented results on total 
cross sections in Ref. and on various differen- 
tial distributions in Ref. ^Sj ■ The latter publica- 
tion also contains technical details of the calcula- 
tion, which is rather complicated. 1 In the follow- 
ing we briefly describe the salient features of the 
calculation and present a selection of numerical 
results that are relevant for a future ILC. 



2. METHOD OF CALCULATION 

The actual calculation builds upon the 
RacoonWW approach 0, where real-photonic 
corrections are based on full matrix elements and 
virtual corrections are treated in DPA. Real and 
virtual corrections are combined cither using two- 



1 Some problems related to the finite width that can appear 
in such a calculation are illustrated in Ref. 1141. 
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cutoff phase-space slicing or employing the dipole 
subtraction method |15| for photon radiation. 

In contrast to the DPA approach, the one-loop 
calculation of an e + e~ — > 4/ process requires the 
evaluation of 5- and 6-point one-loop tensor in- 
tegrals. We calculate the 6-point integrals by 
directly reducing them to six 5-point functions, 
as described in Refs. |16I17) . The 5-point inte- 
grals are reduced to five 4-point functions follow- 
ing the methods of Refs. |17ll8j . Note that this 
reduction of 5- and 6-point integrals to 4-point 
integrals does not involve inverse Gram determi- 
nants composed of external momenta, which nat- 
urally occur in the Passarino-Veltman reduction 
|19| of tensor to scalar integrals. The latter pro- 
cedure leads to serious numerical problems when 
the Gram determinants become small. 

Tensor 4-point and 3-point integrals are re- 
duced to scalar integrals with the Passarino- 
Veltman algorithm as long as no small 
Gram determinant appears in the reduction. If 
small Gram determinants occur, two alternative 
schemes are applied ^21- One method makes use 
of expansions of the tensor coefficients about the 
limit of vanishing Gram determinants and possi- 
bly other kinematical determinants. In this way, 
again all tensor coefficients can be expressed in 
terms of the standard scalar functions. In the 
second, alternative method we evaluate a specific 
tensor coefficient, the integrand of which is loga- 
rithmic in Feynman parametrization, by numer- 
ical integration. Then the remaining coefficients 
as well as the standard scalar integral are alge- 
braically derived from this coefficient. 

As a further complication, also the evaluation 
of the three spinor chains corresponding to the 
three external fcrmion antifermion pairs is non- 
trivial, because the chains are contracted with 
each other and/or with four-momenta in many 
different ways. There are O(10 3 ) different chains 
to calculate. We have worked out algorithms that 
algebraically reduce all these spinor chains to a 
few or, in an alternative method, to a minimal 
set of standard structures without introducing 
numerical problems. These algorithms are de- 
scribed in Ref. [Rj in detail. 

The description of resonances in (standard) 
perturbation theory requires a Dyson summation 



of self-energy insertions in the resonant propa- 
gator in order to introduce the imaginary part 
provided by the finite decay width into the prop- 
agator denominator. This procedure in general 
violates gauge invariance, i.e. destroys Slavnov- 
Taylor or Ward identities and disturbs the cancel- 
lation of gauge-parameter dependences, because 
different perturbative orders are mixed (see, for 
instance, Ref. |2U] and references therein). 

For our calculation we have generalized ^Sj the 
so-called "complex-mass scheme", which was in- 
troduced in Ref. [Jj for lowest-order calculations, 
to the one-loop level. In this approach the W- 
and Z-boson masses are consistently considered 
as complex quantities, defined as the locations of 
the propagator poles in the complex plane. To 
this end, bare real masses are split into complex 
renormalized masses and complex counterterms. 
Since the bare Lagrangian is not changed, dou- 
ble counting does not occur. Perturbative calcu- 
lations can be performed as usual, only param- 
eters and counterterms, in particular the elec- 
troweak mixing angle defined from the ratio of 
the W- and Z-boson masses, become complex. 
Since we only perform an analytic continuation 
of the parameters, all relations that follow from 
gauge invariance, such as Ward identities, remain 
valid. As a consequence the amplitudes are gauge 
independent, and unitarity cancellations are re- 
spected. Moreover, the on-shell renormalization 
scheme can straightforwardly be transferred to 
the complex-mass scheme 

The use of complex gauge-boson masses neces- 
sitates the consistent use of these complex masses 
also in loop integrals. The scalar master integrals 
are evaluated for complex masses using the meth- 
ods and results of Refs. |21I22I2H| . 

In order to prove the reliability of our results, 
we have carried out various checks, as described in 
detail in Ref. [T2]. We have checked the structure 
of the (UV, soft, and collinear) singularities, the 
matching between virtual and real corrections, 
and the gauge independence (by performing the 
calculation in the 't Hooft-Feynman gauge and 
in the background-field gauge 24 ). The most 
convincing check for ourselves is the fact that we 
worked out the whole calculation in two indepen- 
dent ways, resulting in two independent computer 
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codes the results of which are in good agreement. 

3. NUMERICAL RESULTS 

The precisely defined input for the numerical 
results presented in the following can be found in 

Refs. nanj. 

Figure n depicts the total cross section for the 
energy ranges of LEP2 and of the high-energy 
phase of a future ILC, focusing on the semilep- 
tonic final state ud/i~P M . The respective fig- 
ures for the relative corrections 6 to the leptonic 
(shown in Ref. ^21) an d hadronic final states look 
almost identical, up to an offset resulting from the 
missing or additional QCD corrections. Specifi- 
cally, the upper plots show the absolute predic- 
tion for the cross section including the full 0(a) 
corrections and improvements from higher-order 
1SR. The lower plots compare the relative correc- 
tions as obtained from the full 0(a) calculation, 
from an IBA, and from the DPA. The IBA [Hj 
implemented in RacoonWW is based on univer- 
sal corrections only and includes solely the con- 
tributions of the CC03 diagrams. The DPA of 
RacoonWW comprises also non-universal cor- 
rections jS] and goes beyond a pure pole approx- 
imation in two respects. The real-photonic cor- 
rections are based on the full e + e~ — > 4/ + 7 
matrix elements, and the Coulomb singularity is 
included for off-shell W bosons. Further details 
can be found in Ref. 

A comparison between the DPA and the pre- 
dictions based on the full 0(a) corrections re- 
veals differences in the relative corrections S of 
< 0.5% (0.7%) for CM energies ranging from 
s/a ~ f70GeV to 300 GeV (500 GeV). This is 
in agreement with the expected reliability of the 
DPA, as discussed in Refs. [415120) . At higher en- 
ergies, the deviations increase and reach 1—2% at 
y/s = 1— 2TeV. In the threshold region (y/s < 
170 GeV), as expected, the DPA also becomes 
worse w.r.t. the full one-loop calculation, because 
the naive error estimate of (a/ir)(Tw/Mw) times 
some numerical safety factor of 0(1—10) for the 
corrections missing in the DPA has to be re- 
placed by (a/ir)Tw/{y/s — 2Mw) and thus be- 
comes large. In view of that, the DPA is even 
surprisingly good near threshold. For CM ener- 



gies below 170 GeV the LEP2 cross section anal- 
ysis was based on approximations like the shown 
IBA, which follows the full one-loop corrections 
even below the threshold at y/s — 2A% within 
an accuracy of about 2%, as expected in Ref. [S]. 
More results on total cross sections, including 
numbers on leptonic, semileptonic, and hadronic 
final states, can be found in Ref. |12|. 

The distribution in the cosine of the W + pro- 
duction angle is shown in Figure [3for the process 
e + e~ — > ud^~P M at y/s = 500 GeV. Further dis- 
tributions, also for y/s = 200 GeV, are presented 
in Ref. For the W-production-angle distri- 

bution the full 0(a) calculation and the DPA 
agree within ~ 1% for LEP2 energies (see Fig. 12 
of Ref. [Hj), but at 500 GeV the difference of the 
corrections in DPA and the complete 0(a) cor- 
rections rises from —1% to about —2.5% with 
increasing scattering angle (see inset in r.h.s. of 
Figure |2J). Note that such a distortion of the 
shape of the angular distribution can be a signal 
for anomalous triple gauge-boson couplings. 

This work was supported in part by the Swiss 
National Science Foundation. 
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Figure 1 . Absolute cross section a (upper plots) and relative corrections 5 (lower plots) to the total cross 
section without cuts for e + e~ — > udfi'D^ obtained from the IBA, DPA, and the full 0(a) calculation 
(ee4f). All predictions are improved by higher-order ISR. 
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Figure 2. Distribution in the cosine of the W + production angle with respect to the e + beam (r.h.s.) and 
the corresponding corrections (lower row) at y/s = 500 GeV for e + e~ — > ud/i"^. The inset plot shows 
the relative difference A between the full 0(ct) corrections and those in DPA. (Taken from Ref. 
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